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Abstract. - We define and study a three dimensional lattice model which displays a Weyl semi- 
metallic phase. This model consists of coupled layers of quantum (anomalous) Hall insulators. 
The Weyl semi-metallic phase appears between a resulting quantum Hall insulating phase and 
a normal insulating phase. Weyl fermions in this Weyl semi-metal, similar to Dirac fermions in 
graphene, have their lattice pseudo-spin locked to their momenta. We investigate surface states 
and Fermi arcs, and their evolution for different phases, by exactly diagonalizing the lattice model 
as well as by analyzing their topological origins. 
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Introduction. — The interests in topological phases 
of condensed matter have been completely renewed by the 
discovery of the Z 2 topological order, first in two dimen- 
sions (2D) [l}{3] and soon after in three dimensions (3D) 
4 G . This type of order is not described by Landau's 
symmetry breaking theory, but associated with topologi- 



cal properties of electronic structures 7 
bllowing pioneering work in Helium 3 



8 . More recently, 
10 , the focus has 
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shifted to topological order in semi-metallic phases. Such 
a topological semi-metal has been predicted to occur in 
pyrochlore iridates as a result of a combination of strong 
spin-orbit coupling, Coulomb interaction and an antiferro- 
magnetic order preserving inversion symmetry 11 . More 



generally, it has been argued that a semi-metallic phase 
appears in 3D between a normal insulator and a Z 2 topo- 
logical insulator when inversion symmetry is broken 12 
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Proposals following this line include either superlattices 
composed of layers of normal and topological insulators 
subject to time-reversal or inversion symmetry breaking 
[i~3 15 , or magnetically doped bulk Bi 2 Se3 



16] . These 

3D semi-metals, or Weyl semi-metals (WSMs), possess 
gapless chiral surface states and open Fermi arcs termi- 
nating at the projections of the Weyl points [IT , 13 . Such 
unusual features can be manifested through e.g. non-zero 
Hall conductivity 11|13|17 . Recently, realization of nodal 



semi-metals where bands cross along lines has been pro- 
posed [l4] as well as topological semi-metals in fermionic 
optical lattices 
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touch at points (termed Weyl points) where the disper- 
sion relation is linear, in analogy to graphene as a 2D 
(Dirac) semi- metal [19]. The simplest Hamiltonian de- 
scribing the linear dispersion around such a point W is 



the Weyl Hamiltonian %w{q) 



spanned in terms of the 2x2 Pauli matrices er.j, and char- 
acterized by a topological Chern number (that we will 
call helicity by extension of the graphene case) nw — 
sgn [det(Uj_y)] = ±1 20 . The Nielsen-Ninomiya theorem 
imposes that Weyl points must occur in pair(s) with op- 
posite helicity, in any lattice model 21 . When inversion 



symmetry is present, two Weyl points W± located at op- 
posite positions in the momentum space have opposite 
helicity n~w_ = —nw + ■ On the other hand, when time- 
reversal symmetry is present, two Weyl points of opposite 
momenta have the same helicity nw_ = if the a ma- 
trices describe real spin. Therefore, to reconcile the pres- 
ence of both symmetries, two pairs of Weyl points must 
coexist. Such a situation is described by the Dirac Hamil- 
tonian spanned in terms of the 4x47 matrices. The 
conditions for protection of this unstable 3D Dirac semi- 
metal by crystallographic symmetries have been worked 
out in ref. 
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In 3D WSMs, the valence and the conduction bands 



Nevertheless, generic perturbations cou- 
pling the two pairs of coexisting Weyl points open a gap 
- to recover a robust 3D semi-metal, either inversion or 
time-reversal symmetry needs to be broken. 

In this paper, we construct a simple lattice model with 
broken time-reversal symmetry and preserved inversion 
symmetry. This model considers spinless fermions on a 3D 
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Fig. 1: Illustrations of the lattice model with (a) the 3D unit 
cell, (b) the 2D unit cell in the x-y plane, and (c) part of the 
unit cell in the vertical plane which contains both sublattices. 
Nearest-neighbor hopping and next-nearest-neighbor hopping 
are represented respectively by solid and broken lines. Each 
Colored region in the unit cell is threaded by a nonzero mag- 
netic flux: $o /4 out of plane towards the reader for green, and 
$o/2 for purple, with $o = h/e being the magnetic flux quan- 
tum. One specific choice of gauge which is used in this paper 
is also indicated in (b) and (c). 



lattice with properly chosen magnetic flux, and gives rise 
to a topological WSM phase in its phase diagram. Similar 
to the case of graphene, the spinors in our model corre- 
spond to pseudo-spin that originates from the presence of 
two sub-lattices. The Weyl points are obtained without 
resorting to interaction, spin-orbit coupling or magnetic 



order as in the case of real spin 11 . While recent propos- 
als for WSMs were based on perturbing 3D Z 2 topolog- 
ical insulators, the WSM in our model is found between 
a d = 2 + 1 quantum Hall insulator (QHI) and a normal 
insulator (NI) . The d — 2 + 1 QHI was initially proposed 
for layered quantum Hall systems |22p3| , and is obtained 
similarly in our model by layering 2D quantum (anoma- 
lous) Hall insulators. The gapless chiral surface states of 
the d = 2 + 1 QHI lead to a closed Fermi surface that is 
smoothly evolved from the open Fermi arc of the WSM 
through phase transition. We expect this transition to 
manifest itself in the measurement of Hall conductivity. 
In the following of this paper, we start by introducing our 
model and investigating its phase diagram. Owing the 
simplicity of this lattice model, we then examine explic- 
itly the surface bands as well as the resulting Fermi arcs 
in different phases. We provide afterwards a topological 
analysis which accounts for both the bulk and the surface 
properties. Finally we comment on the evolution of the 
Fermi arcs with respect to the global phase diagram. 

Model and Phases. — We consider a tight-binding 
model of spinless fermions defined on a 3D lattice consti- 
tuted of layers of face-centered square lattice with sub- 
lattices labeled by A and B (see Fig. [I]). The nearest 



neighbor hopping between A and B sub-lattices is denoted 
by t, and the second-nearest neighbor hopping between 
A- A or B-B sub-lattice sites is t' . We also introduce an 
on-site energy difference 2A between the two sublattices. 
The planar lattices are coupled with the (nearest neigh- 
bor) inter-layer hopping t± between the A- A or B-B sites 
of adjacent layers. 

The WSM phase arises when certain magnetic flux pat- 
terns are applied to this lattice. For simplicity, we will 
focus in this paper on one of these patterns, indicated in 
Fig. [T] There are two noteworthy properties of this flux 
pattern: it preserves the point group symmetry C^h and 
therefore inversion symmetry, and the flux through each 
surface of the unit cell is zero modulo the flux quantum 
<E>o = h/e. As a consequence of the second property, the 
magnetic Brillouin zone and the original Brillouin zone of 
the lattice overlap exactly, therefore we will not distin- 
guish them in the following. The first property simplifies 
our discussions about surface bands by relating e.g. sur- 
faces parallel to the x-z plane with those parallel to the 
y-z plane. 

With minimal coupling, the effect of the magnetic flux 
is to decorate the hopping terms with additional phase 
factors that are compatible with the flux pattern 24 . This 
permits various gauge choices but different choices do not 
leads to different physics [j] For the specific choice of gauge 
indicated in Fig. [TJd and c, the Bloch Hamiltonian for the 
bulk of the lattice reads 

T~L(k) — (2isinfc ■ a + )a x + (2ts'mk ■ aJ)a y 

+ [A — 2i'(cos k ■ a x + cos k ■ a y ) + 2t± cos k ■ a z ] <r z , (1) 

where a x , a y and a z are the lattice vectors, and a± = 
(a x ± a v )/2. For simplicity we consider a x = a y = a z = 1 
in the rest of this paper. When tj_ 7^ 0, the two 3D bulk 
bands are given by 



E{k) = ± 2 ^(sin 2 fc+ + sin 2 fe_)+ 

t\(mi — vri2 cos k+ cos fc_ + cos k z )' 



1/2 



(2) 



where k± — (k x ± k y )/2, and we have defined two dimen- 
sionless masses mi — A/2t±, m 2 = 2t' /t±. Without loss 
of generality we will focus for now on the range mj > 
and 1TL2 > 0, and discuss the full range at the end of the 
paper. 

The 3D bulk bands are gapped except at points 
where sink + = sinfc_ = and cosk z = — (mi — 
m 2 cos k + cos fc_ ) . The first (two) conditions requires 
{k x ,ky) — (0,0) or (7r,7r). Consequently the last condi- 
tion reduces to 

cosk z — —(mi — m2), for (k x ,k y ) = (0,0); (3a) 
cosk z — -(mi + m 2 ), for (k x , k y ) = (it, tt). (3b) 



^^Note that the inversion symmetry, which can be represented by 
an operator I = <j z for our Hamiltonian JTl , is preserved by an choice 
of gauge, i.e. it satisfies I'H(k)I~ 1 = "Hj—k). 
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Fig. 2: Color-coded 3D bulk band structures of (a) an ideal 
Weyl semi-metal (WSM), (b) and (c) semi-metals when the 
Weyl points merge at G/2 = (0, 0, 0) and (0, 0, n), respectively. 
All the three cases are marked in the phase diagram of the 
model, plotted in (d). The two different Weyl-point-merging 
processes represented by (b) and (c) lead to two topologically 
distinct insulating phases after the bulk gap is opened: one 
features gapless chiral surface states and is a d — 2+1 quantum 
Hall insulator (QHI); the other is a normal insulator (NI). In 
the WSMx2 phase, two pairs of Weyl points coexist. 

Each of the above two equations, when fulfilled with 
double- valued k z e [— it, tt), gives rise to two Weyl points 
related by inversion symmetry and hence having opposite 
helicities (or Chern numbers). 

Let us first suppress the pair at {k x ,k y ) — (tt, 7r) by 
setting mi + TO2 > 1 and validate the other at (fc x , k y ) = 
(0,0) by setting |mi - m 2 | < 1. This is the WSM phase 
regime (see Fig. [2ji) where the bulk bands only touch at 
two distinct points in the momentum space, denoted by 
W± = (0, 0, ± arccos (m 2 — mi)). To examine these two 
points more closely, we expand the bulk Hamiltonian ([I]) 
around them. In the particular case when mi — 7712 = 0, 
the expansions up to linear order in k for W± = (0, 0, ±|) 
and after a global unitary (gauge) transformation H — > 
U^HU with U = exp(— i|u z ) for the purpose of clarity, 
lead to 

H±(q) = v2t q x a x - V2t q y a y =p 2txq z a z , (4) 

where q — k — W± for the two points respectively. These 
are the simplest defining Hamiltonians for Weyl fermions, 
from which one can immediately identify the helicity asso- 
ciated with W+ to be +1 and with W- to be —1 (see Fig. 
[2J1). As |toi — 7772I is increased from 0, the two Weyl points 
move closer to each other in the momentum space, and 
quadratic terms in k z in the Hamiltonian expansions be- 
come more relevant. However, the helicity associated with 



each Weyl point W± remains quantized and unchanged 
due to its topological nature. 

The inversion symmetry forces the merging of a pair of 
Weyl points to occur at G/2 where G is a reciprocal lattice 
vector. This defines the possible boundaries between the 
WSM phase and the bulk insulating phases. Within our 
model the boundaries are given by either mi — 70,2 = — 1 
(the upper boundary) or mi — 777,2 = +1 (the lower bound- 
ary). There, the two bulk bands touch at only one single 
point with vanishing helicity: (0, 0, 0) for the upper bound- 
ary (see Fig. [2]d), and (0, 0, tt) for the lower boundary (see 
Fig.[2j;). Compared with Q, the expansion of the Hamil- 
tonian around the single touch point contains k\a z as the 
leading order term in k z (up to a constant factor), in agree- 
ment with previous studies of merging of Weyl points [12 



or 2D Dirac points [25] . Beyond these transition lines two 
(bulk) insulating phases are reached, denoted by QHI and 
NI in figure|2jl. These two gapped phases are topologically 
distinct, as we will demonstrate below by first examining 
the surface bands and their associated Fermi arcs, and 
then analyzing the topological origin of the surface bands. 

Surface Bands and Fermi Arcs. An advantage 
of a simple tight-binding lattice model is that it facilitates 
explicit investigations of surface bands and Fermi arcs in 
the WSM and its related insulating phases. As a con- 
sequence of the layered structure of our model, a strong 
anisotropy is present : no surface bands exist for surfaces 
perpendicular to the z axis, while for surfaces parallel to 
the z axis we find generically in the WSM phase surface 
bands showing consistent chirality with respect to the z- 
axis. Here we present as examples, in Fig. [3] the surface 
bands for two opposite surfaces parallel to the x-z plane 
(such that k x and k z are good quantum numbers) , for the 
three different phases : the WSM and the two bulk insu- 
lating regimes. As mentioned previously, the results for 
surfaces parallel to the y-z plane can be inferred by using 
the four-fold rotational symmetry that is inherent to our 
model. 

In the WSM case (see Fig. [3^,; we set mi = m 2 = 1 with- 
out losing generality), the surface bands are gapless with 
respect to varying k x for a specific k z g [tt/2, 37r/2]|^]and 
are fully gapped for k z in the complementary range. In 
other words, for a fixed k z in the first range surface states 
sprout from one bulk band (conduction or valence band) 
and merge into the other one as k x is varied; otherwise for 
k z in the second range surface states evolve between two 
parts of the same bulk band. The transition between the 
two types of behavior occurs at the two Weyl points pro- 
jected to the surface Brillouin zone, (k x , k z ) — (0, 7r/2) and 
(k x ,k z ) — (0, 37r/2), where the conduction band and the 
valence band touch each other. Surface states merge into 
the bulk bands exactly at these touch points. An imme- 
diate consequence is the existence of a Fermi arc (for each 
surface band) connecting the two projected Weyl points. 



2 Note that we have shifted the Brillouin zone from 
in the previous discussions for clarity of presentation. 



/hat is used 



p-3 



Pierre Delplacel Jian Lil David Carpentier 2 



m, = 1, m= 1 



w 



+ Fermi 




Fig. 3: Surface bands for three different phases (represented by 
three combinations of mi and m^) corresponding to (a) a Weyl 
semi- metal (WSM), (b) a d = 2 + 1 quantum Hall insulator 
(QHI), and (c) a normal insulator (Nf). Explicitly, we plot 
the two lowest-energy bands for samples built from our lattice 
model that are finite in y but infinite in x and z. The color code 
shows the average position (y) of each state, such that states 
localized at either surface ((?/} ^ y min or y m ax) or extended 
in the bulk ((y) (ymin + ymax)/2) can be clearly identified. 
Gapless chiral surface bands are present in the WSM and the 
QHI, but absent in the NI. The Fermi surface (assuming the 
Fermi energy to be zero) for the WSM appears to be open arcs 
terminating at the two Weyl points, and for the QHI appears 
to be closed lines across the surface Brillouin zone. The three 
scenarios here can be related by examining the evolution of 
surface bands and Fermi arcs upon merging of Weyl points. 



ending the crystal with different sublattices, the shape of 
the Fermi arc may change drastically even when the same 
surface orientation is maintained. 

When the two Weyl points approach each other, the 
band structure deforms continuously, and the Fermi arcs 
of the gapless surface states persist in a certain range of 
the surface Brillouin zone. Depending on whether the two 
Weyl points approach each other through the Brillouin 
zone center or the Brillouin zone boundary, the Fermi arcs 
can either shrink or stretch. It follows that, after the Weyl 
points merge and the bulk gap opens, the Fermi arcs either 



disappear (see Fig. 
louin zone (see Fig. 



z) or persist across the surface Bril- 
3 d) . These two different scenarios for 
the surface states correspond to two different phase transi- 
tions for the bulk (see Fig. [2ji) . The first case corresponds 
to a transition from the WSM to the NI phase, the latter 
being topologically trivial and accommodating no surface 
states. The second case, on the other hand, corresponds 
to a transition from the WSM to the d — 2 + 1 QHI phase 



22 23 where chiral surface states persist while the bulk 
gap opens this QHI phase can be naively expected by 
layering the standard integer quantum (anomalous) Hall 
insulators. Hence in this model, we find the WSM to oc- 
cur between a normal and a topological (quantum Hall) 
insulator. 



Topological Order. — We now turn to an in-depth 
discussion about the origin of these surface states and their 
relation to the Weyl points. Our discussion will be based 
on the dimensional reduction method which involves treat- 
ing one or more momenta, when they are good quantum 
numbers, as parameters of the Hamiltonian.The natural 
starting point for this layered model is to treat k z as an 
effective parameter and reduce the system effectively to 
the 2D x-y plane. This reduction is valid for the bulk 
system or any surface parallel to the z-axis. For conve- 
nience and with a minor abuse of notation, we define a 
fc z -dependent (dimensionless) mass mi(k z ) = m± +cosfc z , 
and rewrite Hamiltonian (JlJ as 



n(k) = g-<r 

with g x — 2t sin k + , g y = 2t sin , 

g 2 = 2t± [mi(k z ) — mi cos k+ cos/c_ 



(5) 



More generally, this observation is valid as long as the sur- 
face is not perpendicular to the z axis such that the pro- 
jected Weyl points do not overlap. The presence of Fermi 
arcs is a robust feature that indeed has a deep topological 
However we emphasize that the specific shape 
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root 

of a Fermi arc is determined by the dispersion relation of 
the attached surface band, which is in turn determined by 
the boundary condition associated with the specific sur- 
face. Therefore by deforming one surface, or simply by 



where er = (cr x , a y , a z ). This 2D model is essentially 
equivalent to the Haldane model for Landau-level-free 
quantum Hall effect [24]. The 2D bulk bands, when fully 
gapped (mi(k z ) =F mi ^ 0), can be characterized by using 
the topological (Chern) invariant defined in the momen- 
tum space as 



c(k z ) = / dk x dk y g • (<9 fcx g x d ky g) , 
47r Jbz 



(6) 
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where g = g/|g|, and the integration runs over the entire 
2D Brillouin zone. For the current model, we obtain 



-U0 -no -no -HO -» 



c(k z 



sgn(M 7r ) - sgn(M )] 



(7) 



with Mq = mi(k z ) — m%, 
M„ = mi(k z ) + m 2 . 



c(k z ) can only take the value ±1 or [24]. For a specific k z , 
a non-zero c characterizes a 2D QHI which accommodates 
one gapless chiral edge mode at its boundary (with the chi- 
rality given by the sign of c) ; while c — signals an NI for 
the effective 2D system with no surface states. The value 
of c(k z ) changes -accordingly a topological phase transi- 
tion happens for the effective 2D model- only when the 
gap closes, that is, when M (k z ) = at (k x ,k y ) — (0,0), 
or/and when M n (k z ) — at (k x ,k y ) — (7r,7r). Hence we 
recover exactly the conditions (|3| for the Weyl points to 
occur in the 3D model. From this point of view, the Weyl 
points appear to be the critical points for the topologi- 
cal phase transitions in a 2D Haldane-like model, with k z 
playing the role of the controlling parameter. It is also 
clear that within a certain range of k z where the reduced 
2D model is topologically nontrivial (cos k z < mi — m\ 
for mi + m2 > 1), edge states persist and form gapless 
surface bands in 3D, while in the complementary range 
of k z (excluding the transition points), the reduced 2D 
model is topologically trivial and a (truly) gapped region 
is expected - the gapless surface bands cease to exist be- 
yond the critical (Weyl) points therefore the Fermi arcs 
end at these spots. Furthermore, as the locations of the 
bulk Weyl points evolve, the range of k z for which the 
(2D) topological order exists can either expand or shrink, 
resulting in eventually the QHI or the NI phase (see Figj3|. 

Next we treat both k x and k y as parameters and reduce 
the system effectively to a ID chain along the z axis. We 
investigate the topological property associated with this 
ID chain in order to account for the absence of nontrivial 
surface bands for surfaces parallel to the x-y plane. In this 
case, the topological invariant is the winding number of g 
-defined the same way as in Eq. [5] for (k x ,k y ) ^ (0,0) 
or (-7T, 7r), but treated only as a function of k z - around 
a given axis when k z runs across the Brillouin zone (i.e. 
from — 7r to 7r). Without entering detailed calculations, 
one immediately finds this winding number to be identi- 
cally zero regardless of the values of k x and k y , by simply 
noticing that g(fc z ) = g(— k z ), hence the trajectories of g 
for k z € [— 7r,0] and for k z 6 [0, 7r] overlaps exactly apart 
from having opposite directions. That is, the effective ID 
chain considered in this case is always topologically triv- 
ial. Therefore no nontrivial boundary states, or nontrivial 
surface bands back to the 3D model, shall be expected for 
this specially chosen orientation. 

Global Phase Diagram. — In the last part of this 
paper, we return to the full phase diagram with respect to 
mi = A/2tj_ and m-2 = 2t' /t±, which has been limited so 
far to the range mi, m,2 > and mi+m2 > 1. We consider 
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Fig. 4: Full phase diagram of our model (the middle panel) 
and two selected trajectories along which the evolution of Fermi 
arcs is examined (the upper and the lower panels). In the phase 
diagram, the grey regions represent Weyl semi-metallic phases 
with only one pair of Weyl points either along (k x , k y ) — (0, 0) 
(hence labeled WSM' ') or along (k x ,k y ) — (7r,7r) (hence la- 
beled WSM M ). Both pairs of Weyl points coexist in the region 
labeled WSM' '"'. The two types of bulk insulating phases, 
labeled QHI and NI, are topologically distinct in terms of the 
Chern numbers defined in the reduced dimensions, given by 
c(k z ) — ±1 and c(k z ) = for every fc z , respectively. Along the 
trajectory (labeled a) passing through the QHI phase, Fermi 
arcs close before reopening, whereas along the trajectory (la- 
beled b) passing through the WSM^ '^' phase, Fermi arcs di- 
vide before rejoining. The helicity of each Weyl point is repre- 
sented by a circle (+1) or a dot (—1). 



several symmetry operations to relate different quadrants 
of the phase diagram. Explicitly, in order to relate two 
points (mi, m.2) and (—mi, 1712) in the phase diagram, it is 
sufficient to transform {k x ,k y ,k z ) — > (k x ,k y ,k z ) + (Tr,Tr,ir) 
and H — > a y T-La y ; in order to map (mi, 7712) to (mi, —1712), 
we use the transformation (k x , k y , k z ) — > (k x , k yi k z ) + 
(ir, 7r, 0) and a x — > —a x ; combined together, the transfor- 
mations {k x , k y , k z ) — > (k x , k y , fc 2 )+(0, 0, it) and a z — > — a z 
map (m 1 ,m 2 ) to (— mi,— m 2 ). 

These symmetry operations immediately extend our 
knowledge of the bulk bands from the first quadrant of 
the phase diagram to the rest (see the middle panel of 
Fig. [ij [^] Namely, the two individual regions given by 
|mi — m-2| < 1 and |mi + 7712I > I are WSM regimes 



3 On the other hand, the boundary conditions, which determine 
the specific details of surface bands and Fermi arcs, are not included 
in this symmetry consideration. 
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with the two Weyl points occurring along (k x , k y ) — (0, 0), 
hence are labeled WSM^ ^; the two individual regions 
given by \m\ + TO2I < 1 and \rri2 — mi\ > 1 are also 
WSM regimes but with the two Weyl points occurring 
along (k x ,k y ) = (71", 7r), hence are labeled WSM' 7r ); the 
semi-metal phase labeled WSM^ 0,71 '' possesses two pairs 
of Weyl points, occurring along (k x ,k y ) = (0,0) and 
(k x ,k y ) = (-7T, 7r); the previously identified QHI and NI 
regimes are now extended by symmetry. 

In this phase diagram, it is instructive to examine the 
evolution of Fermi arcs and their underlying surface bands 
in a coherent manner. This evolution is presented in Fig. [4] 
for a surface parallel to the x-z plane and for two different 
trajectories labeled a and b (see Fig. H). Along trajectory 
a, the open Fermi arc in the WSM^^phase closes before 
opening differently in the WSM^ phase when passing 
through the QHI phase (see the upper panel of Fig. El); 
along trajectory b, the open Fermi arc in the WSM^ 1 ' 
phase divides into two pieces before joining differently in 
the WSMW phase when passing through the WSM^) 
phase [^] (see the lower panel of Fig. [4} . Another notable 
feature in the full phase diagram is the reversing of the 
surface band chirality (with respect to the z axis) and the 
associated sign change of c(k z ) across the mi = line. 
Along the m-i = line, the allowed Fermi arcs appear 
parallel to the k x axis and the underlying surface bands 
are dispersionless in the x direction. Indeed, the sign of 
mi encodes the direction of the magnetic flux penetrating 
through the x-y plane in our original lattice model. 

In this paper, we have proposed a simple tight-binding 
model of spinless fermions which gives rise to a 3D topo- 
logical Weyl semi-metal phase in its phase diagram. We 
have investigated in detail this Weyl semi-metal phase and 
its related insulating phases, especially in terms of their 
topological properties and the resulting surface bands. In 
particular, we have shown that the Weyl semi-metal oc- 
curs between a d = 2 + 1 quantum Hall and a normal 
insulator, in a manner similar to the semi-metal occur- 



ring between a Z2 topological and a normal insulator 12 



We believe that beyond its intrinsic interests, our lattice 
model is a perfect test bench for various properties of the 
Weyl semi-metal phase proposed in various materials, al- 
lowing in particular for numerical studies of its transport 
properties. 

Upon finishing this work we became aware of a preprint 
which focuses on a different lattice model for Weyl 
fermions 126 . 
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